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QUANTUM YANG-MILLS THEORY IN TWO DIMENSIONS: 
EXACT VERSUS PERTURBATIVE 


TIMOTHY NGUYEN 

Abstract. The standard Feynman diagrammatic approach to quantum field theories as¬ 
sumes that perturbation theory approximates the full quantum theory at small coupling 
even when a mathematically rigorous construction of the latter is absent. On the other 
hand, two-dimensional Yang-Mills theory is a rare (if not the only) example of a nonabelian 
(pure) gauge theory whose full quantum theory has a rigorous construction. Indeed, the 
theory can be formulated via a lattice approximation, from which Wilson loop expeca- 
tion values in the continuum limit can be described in terms of heat kernels on the gauge 
group. It is therefore fundamental to investigate how the exact answer for 2D Yang-Mills 
compares with that of the continuum perturbative approach, which a priori are unrelated. 
In this paper, we provide a mathematically rigorous formulation of the perturbative quan¬ 
tization of 2D Yang-Mills, and we consider Wilson loop expectation values on and 
in both holomorphic gauge and in Coulomb gauge with respect to a general metric. We 
establish the equivalence of these two gauges and that both are independent of the choice 
of gauge-fixing metric. Additionally, using holomorphic gauge, we discuss partial results 
showing agreement between the asymptotics of perturbative and lattice Wilson loop ex¬ 
pectations. Our work therefore presents fundamental progress in confirming the paradigm 
that continuum perturbation theory accurately captures the asymptotics of the continuum 
limit of the lattice theory. 
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1. Introduction 

Yang-Mills theory provides a theoretical framework for describing the physics of elemen¬ 
tary particles and has profoundly impacted the study of partial differential equations and 
differential topology. The classical (Euclidean) Yang-Mills action can be written as 

Sym{A) = -^ J {Fa A*Fa) , (1.1) 

where Fa is the gauge-field strength, (•, •) is an ad-invariant inner product on the Lie al¬ 
gebra of the gauge group, e is a coupling constant, and S is the underlying space assumed 
to be a smooth, orientable Riemannian manifold. In the quantum theory, this classical ac¬ 
tion is inserted into a (formal) path integral, from which one can compute various physical 
quantities in terms of a Feynman diagrammatic expansion. The process of evaluating and 
understanding such Feynman diagrams is what leads to many of the basic features of quan¬ 
tum Yang-Mills theory, such as perturbative renormalizability [41, 13, 15] and asymptotic 
freedom [25, 35]. 

On the other hand, quantum Yang-Mills theory can also be formulated on a lattice, 
whereby one obtains a mathematically rigorous, nonperturbative approach that avoids the 
formal aspects of the continuum theory outlined above. Indeed, working on a lattice ensures 
that all quantities can be computed in terms of well-defined finite-dimensional integrals. 
Here, one has to introduce a suitable discretization of the action (1.1) the details of which 
we will return to later. But as one is ultimately interested in a theory that extends all the 
way down to the relevant microscopic scales, one would like to take a continuum limit in 
which the lattice spacing becomes finer and finer. 

Supposing this to be achieved, we obtain two independent, a priori distinct construc¬ 
tions of quantum Yang-Mills theory. While quite different, both the perturbative methods 
of Feynman diagrams and the numerical simulations of lattice methods have yielded spec¬ 
tacular agreement with experimental data in various settings [28, 34, 36]. Naturally then, 
one should consider how these two methods compare at the level of precise mathematics. 
Specifically, since one regards the continuum formulation as perturbative and the lattice 
formulation as nonperturbative, one expects in the limit of small coupling that the two 
formulations should somehow converge. For emphasis, we state this as the following 

Question: For Yang-Mills theory, what is the relationship between the perturbative results 
obtained in the continuum formulation and the nonperturbative results obtained from the 
continuum limit of the lattice formulation, as the coupling constant is sent to zero? 

This paper is an investigation into this basic question, for which we are unaware of any 
prior rigorous work by the mathematical community. 

A priori, our question is well-posed only if we know how to take the continuum limit of the 
lattice formulation. This is a very difficult problem in dimensions three and four, for which 
there exists old work by Balaban [6, 7] that is unfortunately not easily accessible. However, 
we are in the fortunate situation that in two dimensions, quantum Yang-Mills theory has 
a well-known and elegant lattice continuum limit due to Migdal [29] and Witten [43]. The 
aim of this paper is to compare this continuum limit with perturbative two-dimensional 
Yang-Mills theory. 
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Our approach is as follows. The quantities of interest to us are expectation values of 
Wilson loop observables, which are the basic gauge-invariant observables of any gauge the¬ 
ory. Given an oriented closed curve 7 and a conjugation-invariant function / on the gauge 
group G of our theory, we obtain the Wilson loop observable which takes a connection 
A, computes the holonomy of A about 7 , and applies / to this group-valued element: 

Wf.,{A) = f{hol^{A)). 

We can compute the expectation value oiWf^^{A) exactly using the continuum of the lattice 
approach or perturbatively using the methods of Feynman diagrams: 

:= (exact) expectation value 
^fn)pert perturbative expectation value. 

The exact expectation value is defined mathematically precisely in Section 2; for the per¬ 
turbative expectation, defined in Section 3, a few clarifying remarks are in order. We will 
primarily be considering the special case S = 5^. This case is natural for several reasons. 
For the perturbation theory, a compact underlying space is natural since this eliminates 
infrared divergences. Moreover, 5^ is especially simplifying because then there is a unique 
minimal Yang-Mills connection modulo gauge-equivalence, namely the trivial connection^. 
For S of higher genus, the presence of a continuous moduli of flat connections makes the 
perturbation theory more involved. Thus, we take S = 5^ and ^ja)pert denotes the Feyn¬ 
man diagrammatic expansion about the trivial connection. We also consider S = which 
(in the appropriate instances) we can regard as the limit of when the area of the latter 
is sent to infinity. We refer to this limiting procedure as “decompactification”. 

Next, note that in two dimensions, the Hodge star operator * appearing in the integral 
(1.1) is specified entirely in terms of an area form da on T, (and not on a full metric tensor). 
It follows that the coupling constant Aq = has dimensions of inverse area. (In other words, 
scaling the area form by and the coupling constant Aq by preserves the action.) Thus, 
we define the dimensionless coupling constant 

A = Ao|S| (1.2) 

where | • | denotes area with respect to da. Thus, (lTj^.y) is a function of A while (Wf,'i)pert 
is a formal power series in A. 

Finally and quite crucially, the perturbative expectation value (W f,'i)pert I'equires a suit¬ 
able gauge-fixing procedure to be introduced. The most natural choice of gauge to consider 
is Coulomb gauge (also known as Landau gauge). Here, one chooses an auxiliary metric 
g = gij and imposes the gauge-fixing condition d*A = 0to eliminate longitudinal modes. 
The geometric nature of Coulomb gauge makes it applicable for arbitrary S. In the case 
of S'^ however, we have available another choice of “gauge”, namely holomorphic gauge. 
This gauge is referred to as (Euclidean) light-cone gauge in the physics literature. Here, 
writing A as A^dz + A^dz in terms of holomorphic and anti-holomorphic components (with 
respect to a chosen auxiliary conformal structure), holomorphic gauge imposes the condi¬ 
tion = 0. While the interpretation of this condition as a gauge-fixing condition requires 


^We will be focusing only on topologically trivial bundles. This is not a real restriction, see Remark 2.5. 
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some additional analysis [32], Feynman diagrams can be meaningfully generated from this 
ansatz notwithstanding. 

Thus, we consider in this paper 

= (WY>c . 

the perturbative expectation corresponding to Coulomb gauge or holomorphic gauge, re¬ 
spectively, which are defined mathematically precisely in Section 3. Each involves the choice 
of an auxiliary metric compatible with the given area form on 5^ (for holomorphic gauge, 
the data of a conformal structure and the given area form is equivalent to a compatible 
metric). 

Our results can be summarized as follows, which we represent pictorially alongside addi¬ 
tional results from [32, 33] in Figure 1: 

Theorem 1. Consider Yang-Mills theory on {S‘^,da) with arbitrary compact gauge group 
G. Pick any compatible metric for use as a gauge-fixing metric. 

(i) defined using an appropriate regularization scheme, is finite without any 

need for counterterms and is independent of the choice of gauge-fixing metric. 
Moreover, is invariant under area-preserving diffeomorphisms in the sense 

of Corollary 3.17. 

(ii) We have 

(W'/.Jc = (1.3) 

In particular, is also independent of the choice of gauge-fixing metric. 

(hi) Let j be a simple closed curve. Then the series expansion of limA-s.o (Wfa) 
explicitly given by a Caussian integral over q given by the formula (f.7). Moreover, 
this series differs from by exponentially small “instanton” corrections, see 

Remark f.f. 

Conjecture 1. We have 

(1.4) 

Here, ~ means that the left-hand side of (1.4) is described asymptotically by the series on 
the right-hand side^. 

Theorem 2. Let j be a simple closed curve. Then Conjecture 1 holds true up to second 
order in perturbation theory, i.e., up to O(A^), in the decompactification limit S'^ —)■ 

Except for (iii), the above statements and Conjecture 1 readily extend to arbitrary products 
of Wilson loop observables. 

For (i), a regulator is needed since we will obtain divergent Feynman diagrams in the 
expansion defining {Wf^.y)^. We use a heat kernel regulator. That no counterterms are 
needed even as the regulator is removed is shown in Theorem 3.3 via explicit computations. 
Statement (ii), which proves the equivalence of Coloumb and holomorphic gauge (to all 
orders in perturbation theory), is surprising from a purely mathematical point of view 
since the constructions are very different. From a practical point of view, since holomorphic 

^In other words, ifiX) ~ 'l2n>o if i’W — CnA" = o(A^) as A —>■ 0 for every N. We write 

liiriA^o to emphasize that we are considering small A asymptotics. 
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gauge is much more computationally feasible than Coulomb gauge, (ii) represents a dramatic 
simplification. Previously, this equivalence was explicitly checked to second order in special 
cases [8, 9]. Statement (iii) provides an explicit computation of the asymptotics of Wilson 
loop expectation of a simple closed curve to all orders in the coupling. Moreover, it shows 
explicitly how the exact answer and the perturbative answer differ through asymptotically 
zero instanton corrections. There has been previous work on (iii) in the physics literature 
[9, 10, 23], though we were unaware of them when we obtained (iii) independently. The 
result is somewhat mysterious, see Section 5 for some discussion. Theorem 2 provides partial 
confirmation of Conjecture 1, see Theorem 4.1 for a more precise statement. 

Altogether, our results, still to be further completed by resolving Conjecture 1, cor¬ 
roborate the central tenant of quantum field theory which asserts that the Feynman dia¬ 
grammatic perturbative expansion yields an asymptotic series for the exact expectation of 
observables: 

lim (O) ~ (O)perf ■ (1-5) 

Despite enormous efforts to place quantum field theory on firm mathematical foundations, 
establishing the fundamental consistency condition (1.5) in the context of gauge-theories 
appears to have been overlooked by the mathematical community. 



Figure 1 illustrates how our main results above and our subsequent work [32, 33] at present 
fit together. Results (i) and (ii) on S'^ are given by the rightmost arrow and the equality 
between holomorphic and Coulomb gauge in Figure 1. Moreover, when we decompactify 
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5^ to these results about Coulomb gauge and holomorphic gauge carry over to if we 
define the latter gauges as the limit of those obtain from (i.e. take the propagators on 
to be the limit of those obtained on S'^). Additionally, we can consider various axial gauges 
on as we did in [32, 33]. Evaluating perturbative Wilson loop expectations in axial- 
gauge^ provides additional settings in which to test (1.5). One source of ambiguity with 
axial gauge is that it needs a regulator. When one uses stochastic methods (i.e. white-noise 
analysis) to regulate the axial gauge, one obtains “stochastic axial-gauge”, which actually 
decomposes a priori into two cases: partial axial-gauge and complete axial-gauge. In [33], 
we establish the equivalence between these two gauges, hence we have a single stochastic 
axial gauge listed in Figure 1. On the other hand, the Wu-Mandelstam-Liebrandt (WML) 
regulator leads to “generalized axial-gauge” (i.e. the connection has a vanishing component 
along a complexified direction), which is equivalent to holomorphic gauge on [32]. 

The analysis of (1.5) for Wilson loops really splits into three cases; 

(1) decompactify to then compute asymptotics; 

(2) compute asymptotics on 5^ then decompactify; 

(3) compute asymptotics on S'^. 

For case (1), we established (1.5) in the stochastic axial-gauge to all orders in perturbation 
theory [33]. Conjecture 1 therefore concerns cases (2) and (3) using the holomorphic gauge 
(this is denoted by the question in Figure 1). Remarkably, (1) and (2) are inequivalent, 
i.e. Wilson loop asymptotics computed using axial gauge on and holomorphic gauge on 
(the second is the decompactification of holomorphic gauge on S"^) differ, see Theorem 
4.2. A physical explanation for this phenomenon is that the asymptotics on 5^ “remem¬ 
bers instantons” when decompactihed, hence differing from the asymptotics computed after 
decompactification [9]. For case (2), Theorem 2 verifies Conjecture 1 to second order in per¬ 
turbation theory for Wilson loops given by simple closed curves. This result is denoted by 
the O(A^) arrow in Figure 1. Conjecture 1 in cases (2) and (3) can also be verified to all 
orders for highly symmetric configurations [32]. 

Our establishing the fundamental relation (1.5) in stochastic axial-gauge on made use 
of the fact that Yang-Mills theory is (a) free in axial-gauge; (b) has a measure-theoretic 
interpretation in (complete) axial gauge. Thus, relating continuum Yang-Mills theory in 
stochastic axial-gauge to the lattice formulation is to be expected, though by no means 
trivial [17, 33]. (In fact, we know the analysis cannot be trivial because holomorphic gauge 
and stochastic axial-gauge, while both rendering Yang-Mills theory free, yield inequivalent 
Wilson loop expectation values). Thus, what makes Conjecture 1 interesting is that (a’) 
there is no (bona fide) axial gauge on 5^, i.e. a real global direction on which we can 
gauge away; (b’) holomorphic gauge does not have a measure-theoretic interpretation since 
the corresponding propagator does not define a positive-definite pairing. Consequently, re¬ 
lating holomorphic gauge to the asymptotics on the lattice is not at all obvious. Altogether, 
it remains to consider Conjecture 1 for cases (2) and (3) for general loops to all orders. 


^Axial gauge, like holomorphic gauge, while rendering two-dimensional Yang-Mills theory free, leads to 
a “perturbative” Feynman diagram expansion in the sense that the resulting expansion is a priori a formal 
one that does not require the existence of a measure, i.e. an honest expectation. Thus, investigating (1.5) 
is not vacuous in axial gauge [33]. 
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Our paper is organized as follows. In Section 2, we discuss the lattice formulation of 2D 
Yang-Mills theory and describe how its continuum limit yields a rigorous construction of 
a Yang-Mills measure. The most important outcome of this is that one obtains concrete 
formulas for Wilson loop expectation values. In Section 3, we describe the entirely different 
methods of perturbative quantization of continuum Yang-Mills theory. Here we give a rapid, 
self-contained setup of formal perturbation theory in Coulomb gauge using the most direct 
procedure available: the Faddeev-Popov method. While many physics treatments of the 
Faddeev-Popov method use formal arguments to assert that it leads to a gauge-invariant 
construction, proving that this is so mathematically rigorously takes a fair amount of sophis¬ 
tication. For this, we use the Batalin-Vilkovisky method of quantization, which is powerful 
enough to handle the situation in which there are zero ghost modes (which we do find our¬ 
selves in since we quantize about a trivial connection). This is done in Section 3.3, which 
being rather technical, we isolate as an independent section on its own since its importance 
is of a conceptual, not computational nature. We also use the Batalin-Vilkovisky formalism 
to prove the equivalence of Coulomb gauge and holomorphic gauge in Section 3.4. Finally 
in Section 4, we use analytic and Lie theoretic tools to relate the asymptotics of exact 
Wilson loop expectation values to perturbative calculations. We finish with a discussion of 
our results and future directions. Since the methods of perturbative quantum field theory 
are not well-known to most mathematicians, our appendix provides background on Wick’s 
Theorem so as to make this paper as self-contained as possible. 

Acknowledgments. The author would like to thank Vasily Pestun for many helpful dis¬ 
cussions and for pointing out a crucial error in an earlier version of this paper. Pestun 
also referred the author to several references in the physics literature which the author had 
missed. In particular, part (iii) of Theorem 1 appears in a less general form in the previous 
work [10], and it can be generalized straightforwardly to the case of disjoint nested loops 
as in [24, Section 8.3]. The author would also like to acknowledge helpful discussions with 
Greg Moore and with Tom Parker concerning heat kernel methods. 

2. 2D Yang-Mills Measure 

Euclidean quantum Yang-Mills theory, in any dimension, can be given a rigorous formu¬ 
lation on a lattice. As there is no canonical choice for how to discretize a continuum theory, 
many formulations are possible. In two dimensions it is most convenient to work with the 
one due to Witten [43], which extends the original approach of Migdal [29]. This Migdal- 
Witten formulation is invariant with respect to lattice subdivision, so that the continuum 
limit of taking the lattice spacing to zero is in some sense already inherent. 

We recall this formulation (following [26]) in the general setting of when the underlying 
space is a closed, connected surface S endowed with an area form X^da, where da is some 
reference area form and Aq a coupling constant. Let G be a compact Lie group and F 
be a triangulation of S, that is, a finite set of edges (mappings of intervals [0,1] into S, 
injective on the interiors) joined at their vertices such that their complement consists of 
a disjoint union of faces homeomorphic to disks. Assign an arbitrary orientation to each 
edge. To these oriented edges e of F, we assign group valued elements Qe G G, which form 
the basic variables of the discretized Yang-Mills theory. To each face F, we can compute its 
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area |-F| (with respect to da) and assign an orientation to the boundary dF of F, thereby 
inducing an orientation on all the edges which comprise it. Write dF = Fen • • • =t ei as 
a concatenation of edges occurring in their natural cyclic order (well-defined up to cyclic 
permutation), where one has ± according to whether the preassigned orientation of the edge 
Ci is the same or opposite of that induced from dF. Define 

ffOF' ^en * * * ’ 

the product of the group valued elements associated to dF with the appropriate powers. It 
is well-defined up to cyclic permutation of the factors and inversion. 

Pick any bi-invariant metric on G. This is obtained from an ad-invariant metric on its Lie 
algebra g, which we denote by (•,•). One obtains an associated Laplace-Beltrami operator 
A on G. The heat kernel for A is given by the function Kt{g), t > 0, which satisfies 

f){h) = [ Kt{hg~^)f{g)dg (2.1) 

Jg 

for all smooth functions / on G. Here dg denotes normalized Haar measure on G. The 
function Kt{g) satisfies 

Kt{g-^)=Kt{g), Kt{gh) = Kt{hg), for all i/,/i G G. (2.2) 

Assign the weight A'ao|f|(5 '9 f) to F. Properties (2.2) show that this weight is independent 
of the orientation of dF and the cyclic ordering of the edges in dF. The Yang-Mills measure 
associated to P is the measure 

dg.TMda=(y n ■^Ao|f|(5'9f)) dge, (2.3) 

FeF(r) eeE(r) 


on where E(T) and F{T) are the set of edges and faces of P, respectively. The 

Yang-Mills partition function for (S, Aodu) is then 


■^S.Aodo- = / dfir,\odcT- 


(2.4) 


Our lattice action is such that the partition function is independent of the choice of triangu¬ 
lation T. This is a simple consequence of the fact that the heat kernel obeys the convolution 
property 


Etiigig ^)Kt2{gg2)dg = Kt^+t2{gig2) 


IG 


(2.5) 


so that (2.4) is invariant under subdivision. See [43] for further details. 

We can make formula (2.4) more explicit. A surface S of genus h can be represented 
as a 2/i-gon with sides appropriately identified. Applying the above formula using the 
(degenerate) triangulation whose only face is such a 2/i-gon, we obtain 


Zj:,\oda= / ■ ■ ■ahbhaf^^b^^)Y[daiY[dbi. 

JG^h 

In a quantum field theory, we are interested not only in the partition function but also in 
the expectation values of (gauge-invariant) observables. For continuum gauge theories, we 
have the Wilson loop observables Wf^.y{A), which form a rich set of observables since they 
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form a dense collection of functions on the space of connections modulo gauge. This is a 
simple consequence of the following fact; 

Lemma 2.1. Let Ai and A 2 be two eonnections on a prineipal G-bundle P over a connected 
base manifold M. Suppose that their holonomies about every loop based at some given point 
p agree. Then Ai and A 2 are gauge-equivalent. 

Proof. Given any path 7 in M, let Pi{'^) denote parallel transport from 7(0) and 7(1) 
using Ai. Our hypotheses imply that given any path 7 joining p to any other point q G M, 
the automorphism Pi{'^)Pi{'^)~^ of Pq, the fiber of P over q, is independent of 7 . Letting 
q vary, this yields for us a well-defined bundle automorphism. This automorphism is the 
desired gauge transformation intertwining Ai and A 2 . □ 

In the lattice formulation, we would like to compute the expectation (IP/, 7 ) with respect 
to the lattice Yang-Mills measure induced by some triangulation T on S, with IP/,-), suitably 
defined. The lattice formulation makes it clear how to express such an expectation value 
in terms of a closed formula. We assume 7 is regular, namely, that it is a finite union of 
piecewise smoothly embedded curves. Since 7 is regular, we can consider its image as an 
oriented finite graph 7 on S (the maximal components on which 7 is injective constitute the 
edges of T). Embed 7 in some triangulation T of S. Write 7 = ± 6 ^ • • • ± ei as a sequence 
of edges in the order that they occur in the parametrization of 7 , with the ± according to 
whether the given orientation of e* agrees with the one induced from 7 . We then obtain 
9-t = as above. 

Definition 2.2. We have 

(M^/, 7 ) = — [ fi9-y)dtir,\oda- (2.6) 

The subdivision invariance property of our lattice formulation implies that (2.6) does not 
depend on the choice of T containing 7 . 

Strictly speaking, in the lattice formulation of gauge-theories, one chooses a very fine tri¬ 
angulation of S and Wilson loops with 7 adapted to the triangulation. The above formula, 
which allows arbitrary regular 7 (and then extends to arbitrary continuous 7 ) takes into 
account the continuum limit, since as the triangulation gets finer and finer, one can approxi¬ 
mate arbitrary curves. The end result is that (2.6) provides us with an operational definition 
of the Yang-Mills measure. A more refined treatment [26, Definition 2.10.4] shows that such 
a Yang-Mills measure is a measure on P{LT,, G)fG, the quotient space of functions P{L'E, G) 
from LS (the based loops on S) to G, modulo the action of the group Q of gauge trans¬ 
formations. (Given / G F{L'P,G) and g £ G, we have {g ■ f){'j) = 5 '( 7 ( 0 ))/( 7 ) 5 ( 7 ( 0 ))“^.) 
Leaving out many details, such a measure is obtained by constructing a probability space 
in which one has a G-valued random variable P[-f, for every 7 G LS, defined as follows. 
Given 7 (assumed to be regular without loss of generality), the law of is obtained by 
conditioning the Yang-Mills measure on where T D 7 ; namely, one conditions 

on the element g.^^ determined by 7 . Such a law determines a G-valued process 
and thus a measure on G^^ = P{LYi, G), which due to its invariance under Q, descends to 
a measure on F{LT,,G)/G. 



10 


TIMOTHY NGUYEN 


One can informally regard a measure on T{LTi,G) jQ as a measure on the space A/Q 
of connections A modulo gauge transformations. Indeed, we have an inclusion A/Q — )■ 
J-{LTi,G)/Q, given by mapping a connection A to the function hol(.)(A), which takes a 
based loop 7 to hol..y(yl). This makes sense for A sufficiently regular; if A is continuous say, 
then the resulting holonomy function will depend continuously on 7 . Thus, J-{LTi,G)IQ 
represents “generalized connections”. One can also interpret the Yang-Mills measure in 
terms of white-noise measures, with Wilson loop observables being given by stochastic 
parallel transport [17, 37]. 

We only mention these measure-theoretic interpretations so as to note that Yang-Mills 
theory in two dimensions has a rigorous construction in accords with the demands of con¬ 
structive quantum field theory [22]. For our present purposes however, we are only concerned 
with the formula (2.6), which gives the exact expectation value of a Wilson loop observable. 

We will be mainly concerned with the case S = for reasons explained in the introduc¬ 
tion. To that end, let us specialize (2.6) to S = 5^ and the case of a simple closed curve, 
which we will analyze in Section 4. 

Corollary 2.3. Let 7 6 e a simple elosed curve on S‘^, with Ri and R 2 the eonnected eom- 
ponents of S'^ \ ^. Then 

Ig f(9)K^olR^l(g)K^olR2l {g)dg 

^Ao|5^|(l) ^ ^ 

Proof. The graph of 7 consists of a single edge which we label by g. We now apply 
(2.6) with T = 7 and make use of properties (2.2) and (2.5). □ 



Remark 2.4. The above analysis carries over to (equipped with the standard area form). 
For unbounded regions of M^, we use Koo{g) = 1. Using the fact that jQKt{g)dg = 1 for 
all t, the partition function Z]^2^Xod'^x simply becomes unity. Since the area of is not 
normalizable, the effective coupling constant for Yang-Mills theory on is simply A = Aq. 
(We can think of A as Aq times the area of the unit square, which is one.) 

Remark 2.5. The Yang-Mills measure we describe in fact consists of an average of all 
possible topological bundle types over S [27]. In [27], a more refined construction associated 
to a graph G over S and a bundle type a Yang-Mills measure on a configuration space 
that covers G^^^\ However, when A —)■ 0, the dominant contribution to the Yang-Mills 
measure comes from the trivial bundle, since then one has the trivial connection in which 
all holonomies are trivial. (In other words, Kt^gop) is concentrated near ggp = 1 for t 
small.) Nontrivial bundles will be asymptotically zero due to exponentially small factors 
arising from the topology (i.e. instantons). Thus, for our analysis at small coupling, it 
suffices to work with the averaged Yang-Mills measure described above. 

We will compute the asymptotics of (2.7) as A = Ao|5'^| —)■ 0 in Section 4 to give an 
explicit example of the relation between exact asymptotics and perturbation theory. 

3. Perturbation Theory 

In the previous section, we discussed the full quantum Yang-Mills theory, in which one 
obtains expectation values for all possible Wilson loops IF /,7 from a well-defined Yang-Mills 
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measure. This gives a complete, rigorous construction of the quantum Yang-Mills theory 
insofar as it provides a mathematical realization of the formal expression 

f dAe-^YM{A) ’ 

which supposes the existence of a suitable Yang-Mills measure on the space of 

connections (modulo gauge). Here, our basic field H is a connection on the trivial G-bundle 
over S so that the space of connections A can be identified with 17^ (S; g) and the group of 
gauge transformations can be identified with G-valued functions on S. 

In this section, we discuss the perturbative quantization of Yang-Mills theory, whereby one 
computes expectation values not with regard to a true measure but through a perturbative 
expansion in Feynman diagrams about a minimal configuration of the classical action. In 
other words, we proceed by way of the standard paradigm of quantum field theory since its 
earliest inception: regard the right-hand side of (3.1) not as a number but as a notational 
device for generating a formal power series in the coupling constant A. 

The methods by which one generates such a formal power series, and the manner in 
which one establishes its resulting properties, are described with mixed approaches and 
differing degrees of rigor and generality in the physics and mathematics literature. For 
gauge theories, whatever approach one adopts, one ultimately needs to choose a gauge-fixing 
condition and show that the resulting outcome, i.e., the series expansion obtained from (3.1), 
is independent of the choice of gauge. For the case at hand, a degree of sophistication is 
required since we work on curved space, in which case the majority of treatments which 
quantize Yang-Mills theory on flat space do not readily apply. Indeed, flat space techniques 
such as working in momentum space and using dimensional regularization [13, 34] are not 
available to us. 

We thus find it instructive to describe our quantization procedure from first principles, 
albeit in a succinct manner. That being so, this section is written using both rigorous 
mathematics and the physically motivated ideas from which they are derived. We trust that 
the reader finds this two-track narrative enlightening rather than confusing. 

We have two main tasks. The first is to describe mathematically precisely how to gen¬ 
erate the perturbative Feynman diagrammatic expansion of 2D Yang-Mills theory. Here, 
the simplest method is to use the standard Faddeev-Popov procedure. In doing so, one 
introduces a gauge-fixing procedure and corresponding ghosts fields. Moreover, one has 
to choose a suitable regularization and renormalization procedure. On curved space, the 
most natural regularization scheme is obtained from the use of the heat kernels. Using this 
scheme, we will show that in fact, 2D Yang-Mills theory in Coulomb gauge is finite. That is, 
no counterterms are needed to render the perturbative expectation values finite. Physically, 
this result can be anticipated from the fact that 2D Yang-Mills theory is superrenormal- 
izable, so that if one assumes divergences are gauge-invariant (i.e. proportional to |FaP) 
they will in fact be absent by power counting arguments^. We provide a proof of finiteness 
in Theorem 3.3. 

The second task is to show that the formal power series we obtain, which gives a pertur¬ 
bative definition of the right-hand side of (3.1), is independent of the choice of gauge-fixing. 

^Such a naive assumption, however, is unjustified [1, 13]. 
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This second step requires greater sophistication than that involved in the first step, whereby 
we use a blend of ideas from [4, 5, 12, 15] to establish gauge-invariance in Section 3.3. In 
broad strokes, we apply the Batalin Vilkovisky formalism in the form developed by [15], 
which provides a powerful algebraic framework by which to analyze the gauge-dependence 
of a quantization scheme. We find it less enticing to introduce this formalism from the 
start, since for the task of the defining the perturbative expansion and doing computations, 
the Faddeev-Popov procedure is more direct and illuminating. Thus, Section 3.3, while 
conceptually important and necessary, can be treated independently of the others. 

3.1. Faddeev-Popov quantization. To quantize any classical gauge theory, we need to 
choose a suitable gauge-fixing procedure. This is because the path integrals in (3.1) should 
only be over the space of physically distinct configurations, i.e., those which are gauge- 
inequivalent. 

Definition 3.1. A compatible Riemannian metric g = gij is one for which the associ¬ 
ated area form agrees with the given one da on S occurring in the Yang-Mills action 

A compatible metric gij along with the inner product (•, •) on g yields for us an inner 
product fj. {a A */3) on ll*(S;g) and a corresponding adjoint operator d* : n*(S;g) —)■ 
ll*“^(S;g) of the exterior derivative d. We consider the following Coulomb gauge-fixing 
condition d* A = 0, and let 

Ac = {A: d*A = 0}. (3.2) 

Since gauge transformations act via A i—)■ gAg~^ + gdg~^, the above gauge-fixing condition 
eliminates all gauge degrees of freedom (with respect to based gauge transformations) in a 
neighborhood of the zero connection. We call gij a choice of gauge-fixing metric. 

Our path integral 

i J dAWf^^{A)e-^^^^^^ 

is to be replaced with the gauge-fixed path integral 

4 [ dAdet{d*dA)Wf (3.3) 
^ Jac 

The determinant factor det{d*dA) is the Faddeev-Popov determinant that weights gauge 
orbits appropriately^. This determinant can be evaluated via the introduction of anti¬ 
commuting fields, or ghosts. This is because there is a well-defined theory for fermionic 
integration in finite dimensions that produces this determinant factor (see Appendix B), 
and we can extrapolate from this an analogous procedure in the infinite dimensional case. 
We proceed as follows: 

Introduce a pair of ghost fields u},uj, which each are g-valued functions on S. To keep 
them separate, we denote the space of to and Co by n°(S;g) and Q (S;g), respectively. 
Define 

= {uj ■. J ujda = 0} 

®For a rigorous treatment in the finite dimensional setting, see e.g. [31]. 
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and similarly for n_|_. Let 


7^0 


£ = © Q_l © Ac, 


consisting of the total space of gauge-fixed connections and ghosts that are orthogonal 
to constants. The latter condition is to eliminate the kernel of d : n^(5];g) —)■ n^(S;g), 
which arises because of infinitesimal gauge transformations acting trivially on the trivial 
connection. 

We now replace (3.3) with 


where® 


4 f dujdudAWf^^{A)e-^^^'^'^'^ 

A 


<5 = ^ ^ ^ y d*dAU)) da. 


(3.4) 


The integration over uj,uj formally produces the determinant factor det{d*dA) via Lemma 

B.4. 

It is with (3.4) that we can perform a Feynman diagrammatic expansion. This is done 
as follows. Group the extended Yang-Mills action into a quadratic kinetic part and the 
remaining higher order interaction part, which one regards as a perturbation of the former. 
Here, we use Fa = dA + ^[A, A] and dA = d + [A, •]. In this way, we have 

j 


e-5 ^ 


(3.5) 


where 


Skin = ^ J {AA*d*dA) + ^ J {uj,d*duj)da (3.6) 

1 = -^ J {[A, A] A *dA) - J ([yl,yl] A *[yl,yl]) - ^ J {u, d* [A ^]) da (3.7) 


The next step is to write 

4 [ du>dcodAWfAAe~^^^’^’‘^'^ = ^ [ du>du}dAe-^>‘*"WfAA)e^ 

^ Jit ^ Jt 


(3.8) 


and then expand as a formal series in powers of Aq = e^. These terms, multiplied 
with each give multilinear functionals on the space of fields. One then “integrates” 

each of these terms against the “Gaussian measure” ^dujdiodAe~^’‘*'^ defined on €, thereby 
producing a formal power series in A. In reality, what one is really doing is performing the 
algebraic operation given by Wick’s Theorem. This operation is described in Appendix B. 
We describe how this generalizes to the quantum field theoretic setting at hand. 

To invoke the appropriate analog of Lemma B.6, we need to describe the propagator P 
as well as the appropriate expansion of Wf^.ye^ as a Taylor series (i.e. an infinite sum over 
polynomial functions). The Taylor expansion of is automatic from the definition of the 
exponential function and we obtain a formal series in powers of A. For Wf^^, we obtain a 


®It is not necessary to multiply J{oj,d*dAOj) by but this ensures that all terms in the perturbative 
expansion are weighted equally in the coupling constant. 
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Taylor expansion via the representation of hol-),(yl) in terms of path ordered exponentials. 
Namely, 


ho^(A) 



1 + 


where A{t) = Note that this presentation assumes we have chosen an em¬ 

bedding G into the group of unitary matrices U{W) on a vector space W, so that elements 
of g C End(lT) can be multiplied. 

Without loss of generality, we can take / = try to be trace in an irreducible representation 
p : G — 7 - End(E), since the linear span of such functions is dense in the space of conjugation- 
invariant functions on G. This allows us to expand Wf^^{A) as a Taylor series in A: 


WfJA) = trv-(l) + E(-l)" 

n=l 




tTv(^p{A{ti)) ■ ■ ■ p{A{tn))y 


(3.9) 


Here, we also write p : g —)■ End(E) to denote the induced Lie algebra homomorphism. 
The above representation uses the fact that parallel transport is equivariant with respect 
to group homomorphisms: 

p(hol-y(H)) = hol^(p(yl)). 

Altogether, this describes the expansion of the integrand Wf^^{A)e^ as a Taylor series in 
the field variables A,uj,Q. 

Next, the gauge-fixed path integral (3.4) determines for us a propagator P, which is 
the integral kernel of the Green’s operator determined by Skin- Specifically, we have the 
orthogonal decomposition 


^ = im © ker d 
H°(S;g) = ©M 
n°(S;g) =n°(S)©M. 

The decomposition for A depends on our compatible metric, while the other two depend 
only on da. The kinetic action Skin is formed out of the Laplace-Beltrami operator A = 
restricted to £. The Green’s operator we are interested in is the operator 

P = I —0 

' im 

which extends to zero on the orthogonal complement. We have P = P^°^ © pP'^^ according 
to the restriction of P to bosonic A and fermionic uj,u> fields. 

The inner product on H^(S;g) allows us to identify the operator P^°^ with its integral 
kernel 


P^“(x,y) © idg G (im *d © im *d) © (g © g) C (H^(S) © H^(S)) © (g © g). (3.10) 

Here, we have separated the integral kernel of P^°^ as the part P^^^{x,y) which acts as 
A“^ on scalar-valued differential forms and the identity operator on g, which we view as 
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an element of 0 ( 8 > 0 using the inner product on 0 . By abuse of notation, we also denote the 
integral kernel (3.10) simply by It satisfies 

where a, b are 0 -indices with respect to an orthonormal basis of 0 . The fermionic propagator 

pfer 

is the integral kernel of P restricted to © n_|_ suitably interpreted. Namely, we 
have the pairing Ld,uj >—)■ J (©, w) da with which we can use to express the integral kernel of 
A-i : fio ^ 0 ° as an element of © 17^. Since and 17^ are fermionic helds, we pick 
up a minus sign when the order of uj and u) is interchanged. Thus, pl®’’ is defined so as to 
satisfy 

5p/e.cu“(x)a;^(y) = 

This defines pP^ uniquely as an element of A^(f7j_ © 17j^). We have 

p^pbos^pfer 


is the (total) propagator. 

We now suppose = S'^. This way, im = Ac, and A restricted to Ac has no zero 
modes. This allows us to conclude that the Feynman diagrammatic expansion of (3.8), 
following Lemma B. 6 , is formally given by the expression^ 


(W,Cc “ = 


conn^O 


(3.12) 


Here, the subscripts ^‘‘conri’'’ and “0” refer to the fact that we only wish to consider those 
Feynman diagrams which consist of a single component connected® to and which have 
no external edges, respectively. 

The formal definition (3.12) fails a priori because the resulting Feynman integrals we 
obtain are divergent. Thus, we need to choose a suitable regularization procedure, i.e. a 
way of mollifying the integral kernel P to a smooth one P^, e > 0, with P^ —P as e —)■ 0. 
Our regularization procedure is via the heat kernel method, which regulates P via 



im*den“^(s)en° (s) 


(3.13) 


Note that as e —)■ 0, we recover A“^|(r, which is most easily seen by diagonalizing A and 
working on individual eigenspaces. 

Let Pe be the integral kernel of (3.13). It is smooth for all e > 0. We can replace P with P^ 
in (3.12) and obtain a well-dehned formal power series in A. Having chosen a regularization 
procedure as above, we also need to perform renormalization, i.e. counterterms need to be 


there were zero modes, the definition of {WfAp^^rt should be modified to have a residual integration 
over these modes after performing the expansion (3.12). 

®In path integral notation, the normalization factor ^ in (3.8) eliminates disconnected components of 
Feynman diagrams. 
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introduced. These are additional e-dependent (and A dependent) interactions one adds 
to I. One is supposed to arrange the so that 


lim 

e-5>0 


(>Ui 




conn,0 


exists as a formal power series in A. (In general, additional counterterms may also be needed 
to renormalize observables.) A very nice feature of two-dimensional Yang-Mills theory is 
that in fact no counterterms are needed as e —s- 0. This is proven in Theorem 3.3. Thus, we 
can form the following definition: 


Definition 3.2. Fix a compatible metric. The perturbative expectation value of a Wilson 
loop in the Coulomb gauge (3.2) is the formal series in A = Ao|S| defined by 




conn,0 


(3.14) 


This yields for us a mathematically rigorous definition of the perturbative expectation 
value of a Wilson loop, with the preceding discussion revealing its physical origins. In the 
next section, we prove Theorem 3.3, thereby showing that the limit (3.14) exists. We then 
show that (3.14) is independent of the choice of gauge-fixing metric in Theorem 3.4. 


3.2. 2D Yang-Mills is finite. In this section, we prove the following result: 


Theorem 3.3. For a closed surface S, we have 


exists. 


lim e^°^P- {Wf,je^) 


conn,0 


(3.15) 


In two-dimensions, Yang-Mills theory is superrenormalizable, which means that there 
are only finitely many one-particle irreducible® diagrams which are potentially divergent. 
This is a simple inspection: the singularities of P are logarithmic, i.e. P{x,y) ~ log|x — 
y\, so that the type of propagator insertions which lead to divergent integrals is highly 
constrained. Moreover, we only need consider those diagrams that have external edges 
(since these ultimately need to be connected to a Wilson loop observable). 

It is straightforward to see that the only one-particle irreducible diagrams with external 
edges that are potentially divergent as e —)■ 0 are the ones in Figures I-III. Here we have 
omitted powers of Aq appearing in and in /, since for these diagrams they cancel to 

give an overall factor of Ag. Individually, these diagrams are divergent as e —)• 0, but we are 
only interested in sums over Feynman diagrams as occurring in (3.15). What we will show 
is that the sum of these diagrams remains finite as e —)■ 0. The consequence is that the 
sum over all Feynman diagrams in (3.15) become finite as e —)■ 0, since the only possible 
divergences come from a subgraph consisting of the sum of diagrams I-III. 


Observe that each Feynman diagram has two pieces: a Lie-algebraic part and an analytic 
part. This is because the propagator P factors into a differential-form part and the tensor 


Feynman diagram is one-particle irreducible if it cannot be disconnected by cutting an internal edge. 
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Figure I. Tadpole from 4-point vertex: dpb os 



Figure II. Loop from 3-pomt vertex: ^ f ([tI, yl] A 



Figure III. Ghost loop: 

e 


idg G 0 (8> 0. Hence, the integrals arising from Feynman diagrams consist of Lie-algebraic 
contractions and differential-form contractions. For each of the above diagrams, we will 
compute each of these factors separately. (The individual Lie and analytic factors have an 
overall sign that depends on how one expresses Feynman diagrams as Wick contractions, 
but the product of these factors is always well-defined.) Here, a proper understanding of the 
signs and combinatorial factors attached to the different Feynman diagrams is absolutely 
essential, as they affect the sum leading to the cancellation of divergences. Moreover, 
observe that in analyzing the potential divergences of the above diagrams, we can use 

roo 

pbos ^ (3.16) 

Je ni(S;0) 


instead of 
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for the bosonic propagator, where is orthogonal projection onto irnd* using the gauge¬ 

fixing metric g. Indeed, the component of e~^^dt that acts on herd has no affect on 
Figure 2, since the propagator is contracted with the operator d and d^ = 0. For Figure 1, 
the operator Iliind* in is bounded and can be made to act on the two external legs; 
the internal loop is then simply given by That the external legs acquire a projection 

onto im d* is immaterial, since when these legs get contracted into other diagrams using the 
actual propagator P^, the projection imd* has no effect. 

Figure I. (tadpole on 4 point function): 

The only Wick contractions dpbos (—| f ([A, yl] A *[yl, yl])) which are nonzero are those 
which pair yl’s from different Lie bracket terms (since P is proportional to idg). There are 
four such possible contractions. 

Without loss of generality, regard the Wick contraction as contracting the second and 
third copy of A. In what follows, pick an orthonormal basis Ca for g with respect to (•, •). We 
have the Lie algebra factor ([ca, Cc], [cc, e;,]), where Ca and ej, are Lie algebra factors coming 
from the external legs and we sum over repeated indices. Letting C = Ad(ec)Ad(ec), we 
have 

([Ca; Cc]) [Cc) Cfo]) — ([Cc) [Cc; 6a]], e^) 

= Cab 

where Cab is the ab matrix element of C. Thus we have 

(I) Lie = Cab (3.17) 

Note that for G simple. Cab = ~C 2 (Ad( 0 ))( 5 af,, where C 2 (Ad( 0 )) is the value of the quadratic 
Casimir in the adjoint-representation. 

To compute the analytic factor of (I), write the external leg variable as A = A^Ca, where 
A“ is a scalar-valued 1-form, which we can take to belong to im * d. We can treat the 
propagator as a real-valued differential form on S x S, since the Lie algebra factor has been 
accounted for. By abuse of notation, we use the same notation for scalar-valued counterparts 
in this setting. We have 

(/)an = J A (1 (8) *(g) l){P^{x,x) A A^{x)), 

where 4 is the combinatorial factor arising from the number of Wick contractions. We have 
to understand the behavior of P near the diagonal. For x, y G S nearby, we can write 
P{x,y) as a singular part P^^'^^(x,y) plus a remaining smooth part. Since we can work 
with (3.16) instead of P, we consider 

pbos,sing^^^y^ = - (x) (3.18) 

where d{x, y) is the Riemannian distance from x to y (with respect to g) and 

<k^(x)GEnd(T;S,T:S) 

$y(y) = idrjE 
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is smooth in x,y. The metric on S allows us to identify ^y{x) G T*'E<SiT*'E. In particular, 
we can write 

^xix) = Vy,{x) (g) Vy,{x) 

for some auxiliary orthonormal frame v^{x) for T*S. 

For Pg, its restriction to the diagonal is a mollified version of (3.18), namely 


^ (3.19) 

This immediately follows from the fact that ~ for x,y nearby as t —)■ 0. 

In other words, (3.19) shows that the heat kernel regulator (3.16) serves as a short distance 
cutoff, with the 1/2-exponent occuring due to the quadratic relation between time and 
space. So 


{I)a.n — 


1 loge 


- 1/2 


2 27r 

loge“^/^ 
dvr 


/ A°‘{x) ^Vy_{x) ^*{vy{x) ^ A^{x)) 
Jt. 

yl“ A *A^. 


Altogether, the Feynman integral corresponding to Figure I equals 


H) 


{I)Lie{I)c 

loge“^/^ 


dvr 


CabA^ A 


(3.20) 

(3.21) 


Figure II. (one loop graph using 3 point vertices) 


There are two possible types of Wick contractions that lead to divergences. For the first 
type, a propagator joins both copies of dA. For the second, dA belongs to separate Wick 
contractions. The number of possible Wick contractions in each case is d. Call the resulting 
Feynman diagrams (//i) and {II2), respectively. 

For {II2), regard each dA of Ibos as being contracted with the rightmost copy of A in the 
[A, A] in the other copy of Ihos- So we get the Lie-algebraic factor 

{1 12 )Lie — (6c], Crf) ([CbjCdjjCc) 

— ([6<0 6c]) S-d) (6(1, [6fe, 6c]) 

— ([6a) 6c]) [6f,, 6c]) 

= Cab- 

For the analytic factor, it is more notationally convenient to write and then regulate 
afterwards. We have two propagators, which we denote by Pi and P 2 . 


{ihU 




{a^{x) a Pi(x) a *dxP 2 {x)^ {A^{y) A P2{y)) A *dyPi{y)^ 

Sj, 
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To extract the singular part, we restrict to a tubular neighborhood N of the diagonal. Let 
= denote equality up to terms that are finite as e —)■ 0. Using (3.19), we obtain 


(ihU = 


1 1 


4 (27r) 


In L 


yl“(x) A A * 


jx - 

|x - y |2 


A°{y) A A * 


dx'^^{x) A X 

{y - xY^ 


\x - yp 




Only the yi = y 2 terms survive; else we get vi Y ^2 in which case (x — yY{x — yY is odd 
under interchange of x and y. Thus 


{IhYn 


1 1 
4 (27r)2 


(yl“(x) A v^{x)) 


IN 


X - y\‘^ 

X - yY 


{A\y)AvYy)) 


Integrating over y, we see that we get a logarithmic singularity at y = x. When we regulate, 
we obtain 

1 log f h 

(ihr:r = L '' ■ 

For we have a Wick contraction of the two dA's; let the other Wick contraction 

be on the right-most yl’s. We get the Lie-algebraic factor: 


{Ih)Lie = {[ea,ec],ed) {[eb,ec],ed) 
= -Cab 


For the analytic factor, for e > 0, the integral is convergent and one can integrate by parts 
to write f ([A, A] A *dA) = — j {d * [A, A\ /\ A) . This way, we do not have two derivatives 
acting on any of the P^. A similar analysis above shows that the singular part of is 

given by 

1 loff f 

= j aa*a. 

Alternatively, one can also see that the Wick contraction of the two dA’s converges to minus 
the delta function as e —)■ 0, in which case, the integral for {IlYjaY^ equals —{I)an- 
Altogether, the Feynman integral corresponding to Figure II satisfies 




1 log € 

4 ^ 


= 0 . 


CabA'^ A + 


1 log e 
4 




Figure III. (ghost loop) 


Here we have two ghost propagators. One contracts the outermost lo and u, the other 
contracts the inner uj with Q. The Lie-algebraic factor we get is 

{III)Lie = (Cc, [ea,ed]) (crf, [eb,ec]) 

= - ([ea,ec], [eb,ec]) 

= Cab 
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For the analytic factor, we pick up an overall minus sign for the ghost loop. Indeed, 
the ordering of Co and co for the outermost contraction is the reverse of the second inner 
contraction, hence the first propagator picks up a minus sign since it is changes sign under 
interchange of these two fermionic fields. Thus 

= J dP{^''{x) A*A^{x)P^''''{x) j dP^^'^'iy) A*A\y)Pl^''{y). 

Since Pfer{x,y) has singular part \ogd{x,y) near the diagonal, restricting to a tubular 
neighborhood N of the diagonal we get 




2 (27r)2 
1 1 


In 


\x - y\ 




:dx^ A *A‘^(x' 


\x-yr " 'J\x-y\^ 
2 


|x - 2 / 1 ' 


dy'' A *A^{y) 


I l,tz4A“A.Y. 

Jn ‘2 \x- y\^ 


2 (27r)2 

As in case II, when we regulate, we conclude that 

{IIiy:r = log 6-1/2 

Thus, the Feynman integral corresponding to Figure III satishes 

= {iii)ue{iny:r 

= -^ log6-1/2 f c,^A^A*Ay 
471 Js 

Proof of Theorem 3.3: From the above computations, we hnd that 

(I) + {Iiy^^a + {IIiy^'^9 = 0 
Thus, the limit (3.15) is finite since all divergences cancel. □ 


3.3. Metric-independence of Coulomb gauge. In this section, we prove the following 
theorem: 

Theorem 3.4. Let S = 5^. Then is independent of the choice of gauge-fixing 

metric. 

We prove this theorem using the Batalin-Vilkovisky (BV) formalism. The power of the 
BV formalism is that it captures the notion of gauge-invariance in an algebraic manner 
that is well-adapted for perturbative quantization. For additional background and insights 
regarding this formalism, we refer the reader to [15, 38, 42]. We will dive directly into the 
formalism, which is an adaptation of the approach of [15]. 

For our present purposes, a fundamental aspect of the BV formalism consists of being 
able to find a chain complex on which the propagator, viewed as an operator, becomes a 
chain homotopy between the identity map and the projection onto the zero modes of our 
theory. Let us unravel this rather involved statement. Consider the following chain complex 
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degree : — 1 0 1 2 

OO(S, 0 ) ^^"(S, 0 ) (3 22 ) 

X A At xt 

consisting of the ghost field X, gauge field A, antifield At, and antighost Xt. (We use X 
instead of uj for our ghosts now, since in the BV formalism they appear in the action in a 
different form.) The space n^’t(S, 0 ) is just a separate copy of n^(S; 0 ) to keep track of 
the field At. We call this chain complex T, which is our total space of (all) fields, and its 
components have degree listed as above. 

Definition 3.5. A gauge-fixing operator :£—?■£ is an operator of degree —1 that 
satisfies (i) (Q^)^ = 0; (ii) [Q,Q^ is a generalized Laplace-type operator^^. 

A gauge-fixing operator Qt yields for us a Hodge-like decomposition S = im Q © im Qt © 
ker[Q,Qt], More importantly, is allows us to define a Feynman diagrammatic expansion. 
In our situation, we are concerned with the following gauge-fixing operators; 

Definition 3.6. Given a compatible metric g = gij on 5^, define the gauge-fixing operator 
QI via 

f 2 O(S, 0 ) Llht(s, 0 ) 11^(5], 0 ) (3-23) 

Here, Him^ denotes the orthogonal projection onto imd C H^(S, 0 ) with respect to gij. 

Observe that [Q,Q\] = acting in all degrees. On 

% := ker Ag 

is spanned by the constant functions and the fixed area form da on S'^. Thus, Ti is indepen¬ 
dent of the choice of compatible metric g. We can define the pseudoinverse [Q, which 

is zero on ker[Q, Q\\ and the inverse of [Q, Q\\ on the orthogonal complement £±_ to % (the 
complement and the orthogonal projection onto it are also independent of compatible 

a)- 

A gauge-fixing operator allows us to construct a propagator. The propagator we get 
from Qg, which we call the BV propagator, differs from the propagator we obtained in the 
Faddeev-Popov procedure in the previous section. However, BV quantization leads to the 
same set of Feynman integrals as Faddeev-Popov quantization, as we will see shortly. 

Definition 3.7. Given the gauge-fixing operator Q^, we obtain the corresponding BV 
propagator Pg = QI[Q,QI]~^. It is a degree —1 operator on £ which satisfies 

[g,4] = id-n (3.24) 

where H is the orthogonal projection onto P. 

^^This means the operator is of the form V*V -I- F, where V is a covariant derivative, V* its adjoint, 
and F a bundle endomorphism. Also, in what follows we must interpret the commutator [•, •] in the graded 
sense as explained in the appendix. 
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Equation (3.24) is the statement that Pg is a chain homotopy from id and II. The key 
significance of this fact is the following. The right-hand side of (3.24) is ^f-independent. 
Thus, 

[Q, dmetPg] = 0 (3.25) 

where dmet is the exterior derivative on the space of compatible metrics on S'^ (such a space 
is a smooth, connected subvariety inside the space of all metrics) 

Equation (3.25) is the statement that dmetPg is closed as an element Hom]R(T_L), the 
space of M linear maps on £±, endowed with the differential [Q, ■] naturally induced from T. 
On the other hand, since Q is acyclic on T-*-, it follows that [Q, •] is acyclic on Hom]R(T_L). 
Thus, we have 

draetPg ^VCa.[Q,-]- (3.26) 

This is the key identity which allows us to establish gauge-invariance. Observe that Q arises 
from infinitesimal gauge-transformations (and the linearized equations of motion). Thus, 
one can interpret equation (3.26) as stating that the propagator Pg, under changes of the 
metric g, changes by gauge degrees of freedom. Because the underlying classical theory is 
gauge-invariant, this leads to gauge-invariance of the quantum theory. 

The remainder of this section makes the above remarks precise. We need to do the 
following: 

(i) Convert the BV propagator Pg, defined as an operator, into an integral kernel Pg 
so as to be placed on the edges of Feynman diagrams; 

(ii) Describe the BV action so as to obtain the vertices to be used in Feynman diagrams; 
(hi) Use (3.26) and the underlying gauge invariance of classical Yang-Mills theory and 

Wilson loop observables to establish gauge-invariance of Wilson loop expectation 
values. (This exploits the fact that no counterterms, in particular those that might 
have spoiled gauge-invariance, are needed for quantization.) 

We have streamlined our approach in this manner because it then explains what would 
otherwise be many mysterious sign rules in what follows. All such choices of signs can be 
viewed as being carefully crafted so as to ensure (i)-(iii) above hold. In what follows, we 
carry out steps (i)-(iii) in a somewhat abstract, but coordinate-independent manner. Ex¬ 
plicit computations to help make our approach more explicit are carried out in Remark 3.16. 

Step (i); 


To convert an operator to an integral kernel, we need a suitable pairing that induces a 
convolution operator. In other words, given a linear operator K : £ ^ £ and a bilinear 
pairing (•, •)g : T(8)T —)• M, we want to express K as an element^^ K ^ £®£ m. the following 
way. For a G £, we want K{a) to be given by applying I® (•, •)g, with the appropriate signs 
to K®aG£®£®£, which is to say that K is the integral kernel of K with respect to our 
chosen bilinear pairing. 

^^For those hesitant with working on infinite dimensional manifolds, one can always restrict to finite¬ 
dimensional families of metrics. Since all such families are compatible under union and restriction, taking a 
direct limit recovers the full infinite-dimensional family of compatible metrics. 

^^In what follows, all tensor products are completed, see footnote 17. 
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The pairing (•, •)]g and correct sign rule is determined by the following constraint. We 
want 

k{[Q,K]) = Qk{k) 

where k is the map taking an operator to its integral kernel. Here Q acts on S ® S in the 
natural way as a derivation: Q{a 0 b) = Qa 0b + (—( 8 ) Qb. 

Definition 3.8. Define the pairing (•, ■)^y : £ 0 £ ^ by 

(aY)^^, = (aaY) 

where (•, •) is the inner product on g. We then obtain the pairing 

~ ® (3.27) 

Both these pairings are of degree —1, i.e. (a, ^)r / 0 only for |a| + \b\ = 1, and are 
skew-symmetric. 

In particular, (•, ■)^ is an odd symplectic form, i.e., it is skew-symmetric, nondegenerate, 
and of odd degree. From (•, ■)^, we obtain a corresponding convolution operator on £: 

Definition 3.9. Given K £ £ 0 £, dehne the convolution operator K* : £ ^ £ as follows. 
On simple tensors K = Ki0 K 2 £ £ 0 £, we have 

K*a={-l)\^^\Ki{K2,a)^ (3.28) 

This determines * on general K from bilinearity and completion. We say that K is the B V 
integral kernel of K*. 

Note that if \K\ = d, then K* has degree d -|- 1. The sign rules in the definition of the 
BV pairing and in the definition (3.28) are carefully chosen so that the following are true: 

Lemma 3.10. We have that 

(i) Q is skew-adjoint with respect to {■, •)]g, i.e. 

{a,Qb)^ = —(—1)1 I {Qa,b)^ 

(ii) For all K £ £ 0 £, we have [Q,K*] = [QK,*]. 

Proof. These are both straightforward computations. We check (ii): 

[Q,K*]{a) = Q(-l)l^^li^i {K2,a)^ - (-l)l^l+i(-l)""^'^i {K2,Qa)^ 

{QK) * a = (-l)l^^lQi^i {K 2 ,a)^ + {-1)\k^\{-1)\^^\+^K, {QK 2 ,a)^ 

We use the skew-adjointness property (i) to equate the above two lines. □ 
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Given our propagator operator Pg, we obtain the corresponding propagator integral kernel 


Pg £ £ ® £ via 




^9 = ^9* ■ 

(3.29) 

The regulated propagator Pg^^ 

is obtained from regulating [Q,Q\]~^ in Pg\ 



P9,e = Q\ 

(3.30) 

Step (ii): 

Pg,e = Pg,e * ■ 

(3.31) 


In the BV formalism, we have an extended BV action S consisting of the ordinary 
(bosonic) action Sbos, a ghost action Sgh which encodes infinitesimal gauge-symmetries, 
and a Chevalley-Eilenberg action Sce, which accounts for the Lie-algebraic structure of 
the infinitesimal gauge-transformations. These are functions on the space £ in the sense of 
Definition A.l, i.e., they are elements of Sym(T*). Explicitly, 

S = Sbos + Sgh + Sce 
Sbos{A) = Sym{A) 

S,UA\A,X) = j^{A\-dAx)^ 

Sce{X\XuX 2) ^ -\(x\lXuX-2]) . (3.32) 

\ / R 

The quadratic part of S yields a kinetic term 

Skin = ^ I {A,d*dA)-^ I (^A\ dx) (3.33) 

while the negative of the cubic and quartic parts of S yield for us the interaction terms: 

/ = -^/ {[A,A\^*dA)-^ I {[A,A]a*[A,A]) + ^ I (^A^ a[A,X])-Sce. 

(3.34) 

Define Ibos-, Igh-, and Ice to be the first two, the third, and last terms of I, respectively. 

The basis for the extended action is as follows. The BV pairing induces a BV bracket, 
which allows us to convert action functionals to vector fields. Conceptually, it is the (odd) 
Poisson bracket corresponding the BV pairing. 

Definition 3.11. Let E be a local functional, i.e., one given by the integral of a polydif¬ 
ferential function of the fields. Then there is a unique local vector field 5^^F such that 

d22F = {-l)\^\ 

for all V £ £. The BV bracket {F,G} between a local functional F and an arbitrary 
functional G is given by 

{F,G} = dsBVFG. 
where one must interpret d^BYp in the sense of (A.3). 
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The BV action, which has degree zero, satishes the following master equation^^: 

{S,S} = 0. (3.35) 

This equation can be written as 

+ = (3.36) 

since e‘^{Skin, ■} = Q regarded as a derivation on the space Sym(T*) (the action of Q on £* 
itself is the one induced from Q on £ via pullback). In turn, (3.36) decomposes into three 
separate equations, 

~QoIgh Q—l^bos T {^gh} Ibos} — 0 (3.37) 

— Q-lIgh — QiIcE + 2 ^Igh, Igh} + {IcE, Igh} = 0 (3.38) 

{Ice, Ice} = 0. (3.39) 

which depend on one, two, and three ghost fields, respectively. Here Qi corresponds to the 
part of Q which maps £i to £i+i. Equation (3.37) expresses gauge-invariance of the classical 
action Sbos, (3.38) expresses that inhnitesimal gauge transformations act as a Lie algebra, 
and (3.39) expresses the Jacobi identity. 

Thus the master equation encodes symmetries and their algebraic consistency relations 
into a single equation. Expressing all such relations in a compact manner facilitates the 
analysis of symmetries, and in particular gauge-invariance, when quantizing. Next, we de¬ 
scribe the Feynman diagrammatic expansion in the BV formalism, which involves applying 
Wick’s Theorem to the BV propagator and the above BV interaction. It yields the same 
expansion as the Faddeev-Popov prcoedure, as the following lemma shows: 


Lemma 3.12. Fix a compatible metric. Consider the Faddeev-Popov propagator (3.11) 
and interactions (3.7) and the Batalin-Vilkovisky propagator (3.29) and interactions (3.34) 
which for notational clarity we denote by here by , P^^, and , respectively. 

Then 


1- >,odpBV 

hm e e 


,. Ao9pFP 

= hm e e 


As a consequence, we have 
lim 


e^O 




jBV 


= lim e 

conn,0 e^O 


^O^pFP 






jFP 


conn.O 


(3.40) 


(3.41) 


Proof. By direct inspection, the bosonic parts of the propagators and interactions 
coincide. So the only remaining issue consists in comparing the fermionic propagators and 
the interactions that have fermions. We can ignore Ice from the BV interactions because 
the propagator has no component, so that upon setting external leg variables equal to 
zero, those that depend on are annihilated. One can show that the fermionic propagators 
are related as follows. Note that 

pBVM ^ ^ (S; g)) © (S; 0) © g)) 


^^Note that because {•,•} is an odd bracket, it satisfies {E, G} = (—E} for E and G local so 
that the master equation is not vacuous. 
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corresponding to Pg : g) —)■ g) and Pg : g) —)■ g), respectively. On 

the other hand, we know that 




Next, observe that is given by (l(8)*(i) and (*(i(g)l) applied to the two corresponding 

components of This follows, for instance, from the computation 



{lx ® *dy)P^^^^^''{x, y) A A{y) 


P^^'^'^'^{x,y) A*d*A{y) 


= -A-^d*A 
= -P^^A. 


(The implicit g-inner product in the T,y integral has been suppressed from the above no¬ 
tation.) It follows that (1 (8> *d) applied to the 17® Kl 0® component of P^^'Pr eq^^l to 
the n® Kl component of P^^A^r gjg^^ ^.^jg (3_28). On the other hand, if we 

identify *du! G imQs with A^, then 


J {uj,d*[A,uj]) da 


{*dio, [A,(jj\) 



So (the third term of (3.7)) and are equal under this correspondence. One can 
now easily see that (3.40) and (3.41) hold. □ 

Step (hi): 


We now turn to the heart of the proof. We have a few algebraic preliminaries to establish: 

Lemma 3.13. Let v G S. For any linear operator D on £, whieh aets as a derivation on 
Sym(T) and dually as a derivation on Sym(T*), we have 

(i) dov = [dv,D] 

(ii) doK = [Qk, D] 

as operators on Sym(T*). 

Proof, (i) Both dov and [d^^D] are derivations of degree (— l)l^^l, so it suffices to 
check dov = [dv, D] on £*. The statement is then automatic. 
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(ii) Without loss of generality, let K = Ki® K 2 . Using (i), then 

= 2 {^K2dDKi + 

= 2 {dK2[dKi,D] + 

= 2{dK2dKi-,D] 

= [dK,D].U 


Let Kt denote the BV integral kernel of e t > 0. 

Definition 3.14. Define the BV Laplacian to be the “divergence operator” 

A®"" := dKo 

Note that since Kq, being the integral kernel of the identity operator, is a ^-function, 
(F) is not well-dehned for arbitrary F. Moreover, Koo is the BV integral kernel of B. 

Lemma 3.15. If at least one of F or G is a local action functional, then 

{F, G} = A^^{FG) - A^^{F)G - (3.42) 

Proof. This is a straightforward computation by noting that {F,G} is the part of 
dxoiFG) in which the Kq contraction joins F to G. This is what the right-hand side of 
(3.42) expresses, since it subtracts from A^^{FG) those contractions that involve only F 
or G alone. □ 


Remark 3.16. Here, we perform some explicit computations, albeit in a condensed and 
somewhat formal way, to illustrate how our sign conventions work (which we hid behind 
a veil of rather abstract constructions). First, beginning with a hnite-dimensional setting, 
suppose we have a symplectic basis d„.i and of the BV bracket, i.e. 

^ij ~ ) BV ’ 

with even and d^j odd. Then the identity operator has integral kernel 

Kq = (g) d^i + d^i (g) d„.i. 

The divergence operator is then 

A'®^ = dKo = d,,id^i 

and the BV bracket is 


{F, G} = d^iFd,,iG + (-l)l^l9^iF%G. 
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In the Yang-Mills setting, we have the even tangent vectors sx^x) 

tangent vectors , ^x{x) ’ ^ ^ where we have suppressed Lie-algebraic and differential- 
form indices. So 

5 6 5 5 \ 

5A{x) 6A^x) ^ <5Yt(x) 5X{x)) 

and using^^ 



\6A{x)’ SA-^iy) / 
\5X^xy6X{y)/sy 


= 5^‘^\x — y) 
= — y) 


one can check that 


\5Ayx)' 5A{y) / 

\5X{xy5Xyy)/^y 


(3.43) 

(3.44) 


e^{Skin, ■} = Q 

and that the master equation holds: 

(3.37) = -{d* dA, [2l, Y])^ + 

\ / M \ / M 



= 0 


(3.38) = - (^dYi^L(7lt, [., Y2])^+ + - (dYt, [^ 1 ,^ 2 ])^ 

= ((7lt,[d^Yi,Y2])^+o) -(^Ayd[Xi,X 2 ])^-(^Ay[A, [^ 1 ,^ 2 ]])^ 
= 0 

r 

(3.39) = 2[Yi, Y 2 ]—L(Yt, [., Y3])^+ o 

= 2(Yt,[[Yi,Y2],Y3])^+o 
= 0 . 


In the above, <-)• means we symmetrize (in the graded sense) over the odd Y’s and l is the 
contraction of a vector field with a one-form (and the appropriate sign rules must be used). 


^^The minus sign occuring when we switch A and in (3.43) arises from the skew-symmetry of the 
wedge pairing on 1-forms, not from a minus sign occurring in the definition of the BV bracket, as is the case 
for the right-hand side of (3.44). This minus sign later reappears as a minus sign in the third term of the 
identity involving (3.38) worked out below. These minus signs become more transparent if we restore the 
1-form indices we have suppressed. 
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For (3.38), we obtain the vanishing of (3.38) from (Ia. being a derivation with respect to 
[•,•]. The vanishing of (3.39) follows from the Jacobi identity. 

Proof of Theorem 3.4: Let O = he a Wilson loop observable. We want to compute 
the variation of {0)q with respect to variations in the choice of compatible metric g. We 
know that 


^met^g 


Pa = -QPo 


for some Pg ^ £_\_® £_\_ via (3.26) and Lemma 3.13(ii). Since Pg has degree zero and [Q, ■] 
has degree one, Pg has degree —1. However, if we regard Pg as having differential-form 
degree one in the metric directions, then this total degree of Pg is zero. 

It is convenient to combine Pg and Pg into a master propagator 


h = P9 + Pr 

Let =1 denote equality between components of degree one in the metric directions. 
Then we have 

dmet {0)c =1 d^et 

conn,0 

conn,0 

where we are implicitly taking the e —)■ 0 limit of regulated propagators. We have 



In the last line, we used the regulated version of (3.24). Consider the regulated BV diver¬ 
gence operator 

:= YmidK,. 

e-^0 

We have 

(3.45) =1 Ao (-Sgpg -b liin^x, - dK^)Oe^) 

^ connfl 

=1 Ao - dKjOe^) 

conn,0 

=1 Ao ((-go + {/, O}) + + i-QI + Ul, I] + A^^/))e^') 

V / conn,0 

-X^[e^^^P.dK^Oe^) 

conn,0 

We used Lemma 3.13 in the second line and the easily verified algebraic identity A®^(Oe^) = 
({/, 0}-|-A^^0-|-2{L,/}-|-A®^/)e^ in the last line. We have —gO-|-{/, O} = —{Skin^O} = 
0 since O is gauge-invariant. We have = 0 since O has no antifield components. Next, 
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= 0 by the master equation (3.36). We have I = 0 since is symmetric 
in Lie-algebra indices, e > 0, while I is skew-symmetric in them. The final term vanishes by 
the following. The operation contracts da into an entry and 1 into an X entry (and 
there is also the idg component one must contract). We have 8 x^ 0 ^^ — I}oo+dK^I)e^, 

where {•, - joe involves the contraction 9xoo (instead of Oko) so as to connect the input func¬ 
tionals. The term vanishes since iLoo is symmetric in Lie algebra indices. We thus 

have to consider 

-^}oo = {Igh, Ice}oo + -^{ICE, Ice}oo (3.46) 

The second term of (3.46) has an uncontracted argument, for which Pg will not be able 
to contract. Thus the operation |o, which makes external leg variables zero, will annihilate 
all diagrams with external X"^ legs. So we need only consider diagrams arising from the 
first term of (3.46). However, all such diagrams vanish using the condition Pg = 0. (This 

condition follows from (QJ)^ = 0 and the fact that qI commutes with [Q, QI].) Indeed, one 
has an external A and leg in {Igh, Ice}oo, and the placement of propagators Pg on these 
legs implements the operation Pg. Since Pg = 0 for all compatible g, differentiating this 
with respect to the metric is still zero. Thus, when we consider Feynman diagrams arising 
from the total propagator Pg as above, we still get that the contribution from {Igh, Ice}cxi 
vanishes. 

Altogether, we have shown that all terms of dmet {0)(j vanish. This establishes gauge- 
invariance. □ 


As a result of metric-independence, we can deduce that, like the exact expectation, the 
perturbative expectation of Wilson loops are invariant under area preserving diffeomor- 
phisms: 

Corollary 3.17. Let S = 5^. Then for any f and smooth, regular 7 , and ehoice of gauge¬ 
fixing metric, we have 

i/ : 5^ —)• is a diffeomorphism sueh that preserves the area of the connected compo¬ 

nents of 


Proof. Let Cg^da denote Coulomb gauge with respect to a metric g and area form da 
(with g compatible with da). Then {Wf^^o'y)Q ^ can just pull 

back all Feynman diagrams in {Wf^^o'y)^' by the diffeomorphism <I>. Our corollary is proven 
by finding a diffeomorphism 4'^ such that fixes 7 and such that da^ := ^*{^*da) —)■ da 
uniformly (in with /c > 2) as e —)■ 0. This is because by choosing a gauge-fixing metric 
ge compatible with da^ such that g^ ^ g uniformly, we have on one hand 


fAoi) Cg^da ~ 


9*('S>*{g)),dae 




9e,d(Te 


where we used independence of the gauge-hxing metric when the area form is fixed, and on 
the other hand 




Cg^ ,d(7 g 


{Wf, 


■'y/c, 


9,do- 


0 


(3.47) 
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as a formal series in A. Indeed, the propagator evaluated at distinct points depends con¬ 
tinuously on the gauge-fixing metric, so that the renormalized integrals occurring in (3.47) 
depend continuously on the gauge-fixing metric and area form. 

To construct we proceed as follows. Thicken the image of 7 to a “ribbon” of width 
e (with respect to some arbitrary fixed metric), whose boundary consists of smooth disjoint 
simple closed curves of distance e from 7 . Let be a d-thickening of dR^, with 6 < e, so 
that it is a union of annuli, one for each component of dRe- 

We can define a diffeomorphism T : 5^ —)• 5*^ such that (i) it maps R^ and the compo¬ 
nents of 5^ \ Re to themselves; (ii) is the identity outside of Ns; (iii) = da on N^s 

for c < 1 small (simply choose det T appropriately on Nf.s)- Consequently, each component 
D of 5^ \ dRe is such that its area with respect to da and T*((icj') differ by 0 (e) and 
da = ^*(da') near the boundary. By the work of Moser [30], on D we can homotope the 
diffeomorphism T : D —)■ D (with the homotopy being identically one near dD) to Ti such 
that 'I'][(dcj') = da' + da(e), where da(e) is an area form supported away from the boundary, 
is of order 0(e) in C^, and is such that the total area of D with respect to "^Kda') and 
da' + da{e) agree. Doing this on each component of 5^ \ dR^, we have constructed the 
desired □ 


3.4. Coulomb gauge = Holomorphic gauge. In this section, we introduce the holomor- 
phic gauge and use the Batalin-Vilkovisky formalism of the previous section to prove its 
equivalence with Coulomb gauge. The interpretation of holomorphic gauge as a gauge-fixing 
condition is subtle, see [32]. Nevertheless, holomorphic gauge is defined as follows. 

Pick any conformal structure on S'^. (This, together with the area form is equivalent 
to a choice of compatible metric.) Using the resulting complex structure it induces, we 
can complexify the space of connections A = D^(S;g) to Ac = D^(S; 0 c) where 0 c is 
the complexification of 0 . The Yang-Mills action extends complex-linearly to connections 
belonging to Ac (by extending the inner product on g complex-linearly). We say A is in 
holomorphic gauge if it is a differential form of type ( 1 , 0 ), i.e. it belongs to D^’°(S;gc). 
The Yang-Mills action in holomorphic gauge becomes 

YM{A) = j{AAd*dA), A e gc). (3.48) 

Indeed, the quadratic terms in the curvature Fa vanish in holomorphic gauge. 

On since there are no nontrivial holomorphic 1-forms, the pairing (3.48) is nonde¬ 
generate. Hence, the operator B * B : D^’®(S) —)■ D®’^(S) is invertible and it has an integral 
kernel, with respect to the wedge pairing, belonging to D^’*^(B) KlD^’®(S). This yields yields 
a corresponding holomorphic gauge propagator 


by tensoring with the identity tensor in g ( 8 > 0 . Explicitly, if z and w are local holomorphic 
coordinates on C = 5^ \ { 0 } with respect to the standard conformal structure on 5^, we 
have 


Phoi(z,w) 



z — w 


z — w 


dw 




(3.49) 
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From the holomorphic gauge propagator, we obtain^® 

(3-50) 

0 

We have the following result; 

Theorem 3.18. Pick any compatible metric on S‘^ determining a corresponding Coulomb 
gauge and holomorphic gauge. Then 

{lU.-,)o = (WY>»<,r (3-51) 

Proof. In the same way that holomorphic gauge needs to be interpreted in terms of 
a real integration cycle [32], the proof of Theorem 3.18 exploits this idea in the Batalin- 
Vilkovisky context. Namely, we complexify the BV complex (3.22) and connect the Coulomb 
gauge and holomorphic gauge through a one-paramemter family of gauge-fixing operators, 
one which interpolates between the subspace im * d c and a totally real subspace of 
Oho, 

We have the decompositions 

= imd © im * d 

with and the ©i eigenspaces of *, respectively. Concretely, we have 

nho = i{-dg + *df)} 

given by the graphs of ±i* : —)■ where / and g denote real-valued functions. Thus, 

we can define a totally real subspace of by restricting the graph of i* to im * d; 

;= {a G : a = *df — id/}. 

We have 

flho ^ j^l,0 ^ -^1,0^ 

Next, we have the following complexification of the BV-complex: 

(3.52) 

The differential Q and the Coulomb gauge-fixing operator Q\ defined as in (3.23), extend 
complex-linearly. Denote the complex (3.52) by Sc, it consists of the terms S^ supported in 
degree k, —1 < k < 2. 

Define the one-parameter family of chain isomorphisms Ut '■ Sc ^ Sc, 0 < t < 1, hy 

JJ \ _ f riiind + (1 + fi*)nim*rf A: = 0, 1 

■ \id k^l 

where the projections Ilimrf and Ilim^d are complementary orthogonal projections on 
(with respect to the chosen compatible metric) extended complex linearly. Indeed, one can 

^^Note that because there are no interactions, all Feynman diagrams are automatically connected. 
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check that Ut commutes with the differential Q. We have Uo is the identity and U := Ui 
maps im *d C to 0^’°, k = 1,2. 

Define 

q\ = UtQ^Ui\ 

Then 

[Q,Ql] = Ut[Q,Q^Ur^ = [Q,Q^] 

since A = [Q,Q^ commutes with Ut- We have the t-dependent propagator 

Pt = Ql[Q,Q^-^ = UtQ^[Q,Q^]-%-\ 

It satisfies 

[g,Pi]=id-n (3.53) 

just as we had in (3.24) with the Coulomb propagator. 

Hence, the exact same proof as the proof for Theorem 3.4, with the e-regulator inserted 
and then taken to zero, shows that Wilson loop expectations with respect to the propagators 
Pt (which are the integral kernels of Pt) are independent of t. Since Pq is the Coulomb gauge 
propagator, it remains to check that Pi is the holomorphic gauge propagator. Indeed, if we 
do this, then 

e^odp^Oe^ =e^°^PhoiO . 

conn,0 0 

for a gauge-invariant observable O, since I vanishes whenever all the entries belong to 
Moreover, for Pi, no regulator is in fact needed since the holomorphic gauge propagator 
(3.49) is bounded as we approach the diagonal, which renders (H7,7)/jo/ finite. 

So we need to check that Pi inverts B * B : —)■ gg a = Bf ^ for some 

arbitrary function /. We have 

UQ^U-^[Q,Q^]-^B*Ba = UQ^U-^[Q,Q^]-^B*BBf 

= -^UQ^u-^Bf 

= -\uQ^u-^ Q(l-«)#) 

= \u{*df) 

= ^(1 + ^*)(* d /) 

= i* Bf 
= a. □ 
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4. Exact Asymptotics vs. Perturbation Theory 

The relation (1.5) for A(/>^ theories in dimensions two and three have been established, 
where it is known that n-point correlation functions have asymptotic expansions in small A 
equal to the formal series one obtains from perturbation theory (in fact, the latter can be 
Borel resummed to recover the n-point functions exactly) [22, Ch 23.2]. Such results do not 
directly address our line of inquiry however, since for scalar field theories, both perturbative 
and nonperturbative calculations involve the same scalar field and therefore use identical 
formulations. What makes the investigation of this paper notable is that it compares two 
different formulations: lattice (group-valued fields) versus continuum (Lie-algebra valued) 
fields. 

For Yang-Mills theory, the central difficulty is that while lim;v_^o ™ two-dimensions 

can be evaluated by determining the asymptotics of heat kernels on Lie groups, the compu¬ 
tation of (fl^/, 7 )pgrt) to all orders in A is comparatively much harder to perform. Indeed, for 
Coulomb gauge, as the order in A increases, one has an increasingly complicated set of Feyn¬ 
man integrals to calculate. While holomorphic gauge is simpler since the theory becomes 
free (i.e. the interactions I do not contribute to Feynman diagrams), the integrals one has 
to contribute are still highly nontrivial. This arises from G being nonabelian (the non¬ 
trivial case), since then the combinatorics of Wick contractions arising from Lie-algebraic 
insertions into Wilson loop operators complicates the analysis [32]. 

Nevertheless, we are able to obtain the following results. First, let us make more precise 
Conjecture 1 to include the decompactification limit. We use subscripts 5^ and to denote 
the space on which exact or perturbative expectations are computed^®. 

Conjecture 1’. We have 

(4,1) 

and 

Ao“o ~ iWf,^)R^hol ■ (4-2) 

Here hm 52 _^]R 2 denotes the decompactification limit in which the area form on the sphere 
increases and limits to the standard area form on M^. In doing so, we regard 'y as a fixed 
subset of S‘^\{oo} = In (4-2), the right-hand side is by definition a formal power series 
in Aq while the left-hand side is the limit of the Aq —)■ 0 asymptotic series of g 2 as 

Thus, (4.1) and (4.2) correspond to cases (3) and (2), respectively, from the introduction. 
We now prove Theorem 2 which says that up to second order, the decompactification of the 
asymptotics on 5^ can be captured by holomorphic gauge on M^, thereby providing partial 
confirmation of (4.2). 

Theorem 4.1. Let j be a simple closed curve. Then 

d'i ~ + 0{4). (4,3) 

^®On (equipped with the standard area form), our dimensionless coupling constant A equals Aq. We 
can regard A = Ao|[0,1] x [0,1]|, where [0,1] x [0,1] is the unit square. 
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Proof. In [32], we computed the right-hand side of (4.3) for 7 a simple closed curve. 
The result was given by a simple Gaussian matrix integral. On the other hand, Theorem 
4.3 shows that the left-hand side of (4.3), to all orders in Aq, is also given by the same 
Gaussian matrix integral. □ 

The next result indicates a subtlely in the various limits and gauges that arise in two- 
dimensional Yang-Mills theory; 

Theorem 4.2. We have 

lim lim {Wf ^ hm (ITf ) 2 , (4.4) 

i.e. decompactification and small coupling asymptotics do not commute. Consequently, 
stochastic axial-gauge and holomorphic gauge are inequivalent on . 

Proof. Let 7 be a simple closed curve. The left-hand side (4.4) is computed to all 
orders in (4.8). The right-hand side is determined from 

= [ K,^^R^ig)f{g)dg 
Jg 

= (e-"ol^l^/2/) (1) 

^ g-Ao|R|c 2 (p)/ 2 ^-^^ ( 4 . 5 ) 

where / = trp and C 2 (p) is the quadratic Casimir for the irreducible representation p. The 
functions (4.5) and (4.8) are not equal. In fact, they define entire functions in Aq, and so 
their series expansions about Aq = 0 are not equal. Example 1 at the end of this section 
shows that (4.5) and (4.8) can disagree at second order in Aq. Hence by Theorem 4.1, on 
stochastic axial-gauge (which computes the exact expectation (IT/,-y)g 2 via [33]) and 
holomorphic gauge are inequivalent. □ 


We conclude with some explicit computations of the full small coupling asymptotics of 
(Wf,'y)s^ for 7 a simple closed curve. The final answer is remarkably simple. 


Theorem 4.3. Let 'y be a simple closed curve on 5^. Define 


P ^ |^ 2|2 


(4.6) 


where Ri and R 2 are the two regions in the complement of 'j. Then for any compact Lie 
group G, we have as A —)■ 0 


lim ~ 1 ^ /(exp(X))e-l^l72->AY, (4,7) 

In the decompactification limit (in which ji? 2 |/|<S'^| —^ 1; where R 2 is the unbounded region 
of'y viewed as a subset we obtain 


lim lim (Wf-,) 

5'2^]K2 Ao-S>0 


52 






(4.8) 
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Proof. Let iL be a maximal torus of G with Lie algebra f). For the time being, we 
normalize our inner product (•, •) on g so that the volume form it induces on g is the 
normalized Haar measure on G. This is so that the normalized Haar measure dg occurring 
in the definition (2.1) of the heat kernel Kt coincides with the Riemannian volume induced 
from the bi-invariant metric determined by (•, •). (Otherwise, we will pick up awkward scale 
factors in what follows.) Let dY and dX denote the volume forms on g and f) induced from 
the Haar measure on G and H, respectively. 

Let L be a regular element of H. We have the exact formula [11, Section 7.2] 


Kt{h) 


E 

YSexp“l(h) 


1 

(277^)*^/^ 


(dety(expj) 


1/2 



st 

12 


(4.9) 


which expresses the heat kernel on a compact Lie group as a sum over geodesics. Here, s is 
the scalar curvature of G, defy (exp^) is the determinant of the differential of the exponential 
map exp : g —)■ G at T G f), and | • | is the norm with respect to the inner product on g. 
This determinant can be written as 

<ie.Uexp.) = 4^, (4.10) 

where j and J are ad-invariant functions on G and g given by 

i(exp(y))= l[ |e“(^)/2 _ e-“(^)/2|2 
oeR+ 

J{Y)= n \a{Y)\\ TGf). 
oei7+ 


Here, R+ is the set of positive roots of G. Thus the leading t —)■ 0 asymptotics of Kt{h) for h 
near 1 is given by the term of (4.9) with Y of smallest length, in which case [Tj = dist(/i., 1). 

Let ti = so that -|- ^2 = A. Let W denote the Weyl group of G. Then by the 

Weyl integration formula 

fid) Kt,{g)Kt^{g)dg = J^f{h)Kt^{h)Kt2{h)j{h)dh (4.11) 

where dh is Haar measure on H. The asymptotics of (4.11) is determined entirely by the 
integral in a neighborhood of 1. Using (4.9) and keeping only the leading shortest geodesic 
term, (4.11) as A —)■ 0 is asymptotically equal to 


|lT|(27r)'^/2tf^t2/^ 


_izP _llP 
/(exp(y))e 2*1 e 2*2 


( j(exp(y)) 

V J(Y) 


j{ex.p{Y))dY 


1 

|lT|(27rp)‘='/2 


/(exp(y))e-l^l"/2pj(y)^y^ 


On the other hand, this last integral is just 

by the Lie algebra version of the Weyl integration formula. 


(4.12) 


(4.13) 
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Suppose now our inner product on g is not the Haar inner product but (? times it. Ob¬ 
serve that this scaling can be absorbed into the * appearing in the Yang-Mills action, which 
has the effect of scaling areas by c“^. In particular, we replace p ^ pjc? in (4.13). This 
is equivalent to replacing the Haar inner product with (? times it, i.e., our chosen inner 
product. Thus, (4.13) holds for the general case. □ 


Remark 4.4. In the proof of (4.7), we neglected exponentially small contributions to 
(4.9) arising from non-minimal length geodesics joining the identity element of H to h. 
Such geodesics have nontrivial winding and we can regard their contribution to Kt{h) as 
“instanton contributions”. 

In [23], additional computations are done that provide an explicit check of Conjecture 1’ 
to all orders in A, namely the case of products of Wilson loops formed out of concentrically 
nested circles. The author has also done nontrivial partial checks in the case of a figure 
eight loop to second order in A. 


Example 1. Let G = SU{2) equipped with the bi-invariant metric induced from trace in 
the standard representation on g. We want to compute the exact asymptotics of 
on 5^ where Xm is trace in the m-dimensional representation. Define the function 

F{p) = e-P/\2-p). (4.14) 


As A —)■ 0, we have 






F((m - l)^p) + F[{m - S^p) H-h F{p) 

Fl{m — l)^/o) -I- F((m — 3)'^p) F{2^p) + 1 


m even 
m odd. 


(4.15) 


Indeed, we use (4.13) or (4.12) with 

/ = Xm(exp(0/)) = + ... + 


where we have the generators 



forsu( 2 ). So (4.12) becomes 

[ Xm(exp(A))e“l^l'^/^^(iA = f Xm{ex.p{xil + X 2 J + ^ ^ , 

Jsu{2) 4R3 [TTpy/^ 

A 

= X™(exp(W))e-l^l'/^02^F 

This yields (4.15). 

On the other hand, on we have {Wf^^) is given by 

/ Kxo\R\{g)Xm{9)dg = e 4 rn (4.16) 

Jg 

since Axm = Am- In particular, the 5^ —)■ M? limit of (4.15), in which p becomes 
AojFj, does not equal (4.16), beginning at second order in Aq. 
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5. Discussion and Further Directions 

We conclude with some natural questions and directions for future research. 

1. Why do we obtain asymptotic series that are entire functions? The explicit asymptotic 
series we computed in Section 4 all defined entire power series (for / a polynomial function 
such as trace) in the coupling constant. Yang-Mills theory being free in axial and holomor- 
phic gauges may make this consequence seem natural. On the other hand, the philosophy 
of resurgence theory [18] predicts that our series expansions should be non Borel-summable 
owing to the presence of instantons. At present, we have no way of reconciling this predic¬ 
tion with what is actually the case. 

2. Compute asymptotics of more complicated Wilson loops. The asymptotic formula (4.7) 
and its generalization to products of Wilson loops obtained from nested, non-intersecting 
simple closed curves [24] suggests that perhaps similar asymptotic formulae can be obtained 
for Wilson loops involving more complicated curves. What is remarkable about these for¬ 
mulas is that they are given by simple Gaussian integrals over the Lie algebra that are of 
the type appearing in random matrix theory. It would be of interest to investigate how 
complicated one may make a Wilson loop and still continue to obtain such formulas for 
exact asymptotics (to all orders). 

3. Consider more general topologies. In extending our results on the independence of the 
choice of gauge-fixing to surfaces of higher genus, the Batalin-Vilkovisky formalism would 
have to be carried out in the case when one has zero modes for the bosonic kinetic operator 
(corresponding to having a continuous moduli of flat connections). Since this difficulty is 
a finite-dimensional complication orthogonal to the infinite-dimensional nature of quantum 
field theory, we do not anticipate any fundamental obstacles. Perhaps a very interesting 
scenario to consider would be to consider compact surfaces with boundary. Here, boundary 
conditions come into play and we have not considered how the two formulations, exact and 
perturbative, line up in this regard. 

On the other hand, in higher genus, relating the asymptotics of the exact expection to 
perturbation theory may prove to be difficult, since now the set of flat connections about 
which to do perturbation theory becomes nontrivial. For each flat connection, we obtain a 
corresponding propagator formed out of the Green’s operators for the Laplacian twisted by 
the such a flat connection. It is unclear to what extent explicit computations can be done 
for nontrivial flat connections. 

4 . Find a priori analytic relationships between the lattice and continuum formulation of 
quantum gauge theories. This is of course not a new question but a very difficult one, since 
the basic variables in the two formulations are of different natures (Lie group elements versus 
Lie-algebra valued 1-forms). So while an infinitesimal group element being approximated 
by an element in the Lie algebra provides the basis for how one discretizes the continuum 
theory, without precise estimates, one cannot establish any clear relation between the two. 
A resolution of this question would presumably shed light on some of the above questions 
and observations we have made. 
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Appendix A. Graded Vector Spaces 

A (real) graded vector space V is a vector space together with a decomposition V = 
into vector spaces V) in degree i. If n G V, then |r;| = i denotes its degree. An element is 
even or odd according to whether its degree is even or odd, so that we have a corresponding 
decomposition of V = Vev © Vodd by parity. An ordinary vector space yields a graded vector 
space concentrated in degree zero. 

For ordinary vector spaces V, one has the familiar notion of Sym"'(I/) and A"'(V), the 
symmetric and exterior powers of V. For graded vector spaces, one defines symmetric 
powers in the graded sense. Namely, let ©"'V be the graded vector space whose graded 
components are 

(©-V), = 0 V, ©•••©v„. 

*lH- \-in=i 

We have an action of Sym„ such that a transposition of adjacent elements acts via 

Then Sym'^(V) is the Sym„-invariant subspace of ©""V with respect to the above action. 
We write 

Sym(V) = 0Sym”(I/) 

n>0 

to denote the total symmetric algebra on V. Thus Sym(V) is a symmetric algebra and an 
exterior algebra (in the ordinary sense) on Vev and Vodd-, respectively. 

If V is a graded vector space, then its dual space V* is the graded vector space given by 

vr = {v-^r, 

that is, the degree i component of V* is the dual space of V-i. In this way, the evaluation 
pairing 

V* © V ^ M 

is a degree zero map. 

Definition A.l. Given a graded vector space V, a function on V is an element of Sym(V*). 
A vector field on V is an element of Sym(I/*) © V. 

For V an ordinary vector space, the above coincides with the ordinary notion of a poly¬ 
nomial function or polynomial vector field. 

A (linear) map / : V —)■ W of graded vector spaces has degree p if /(Vi) C Vf_^p. If the 
degree is not specified, it is understood to be degree zero. The commutator of two maps 
f,g :V ^ V of degrees |/| and \g\ is defined using the appropriate sign rule: 

[f,g] = fg- (-i)'^"^'^/- 

A.l. Directional derivatives. The space (f/*)®^ denote the space of n-multilinear maps 
from V^^ to M. It has a natural action of Sym„ induced from the one on For v G V, 

define the contraction operator 

Ul © • • • © U* l-A ^(— + + (g) . . . (g) y*__^ (g) y*^^ (g) . . . (g) y^ _ (A.l) 
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In other words, dy is the directional derivative with respect to v, where in the graded setting, 
it is a derivation of degree |n| using the above signed Leibniz rule. 

More generally, given an element K = u 0 v G ^ we can define the operation 

Qk = ]^dydu. (A.2) 

This operation extends bilinearly to any K G and depends only on the component of 
K in Sym 2 (’L). 

One can consider contractions using “non-constant coefficient” vector fields, i.e., elements 
of Sym(y*) (g) y. If fv G Sym(F*) (g) y, with / G Sym(F*), then 

dfv = fdv 

Likewise, we have 

dyf = (A.3) 

using the usual sign rules. 

A.2. Infinite-dimensional case. The above considerations generalize to the infinite-dimensional 
setting needed for quantum field theory. Instead of a finite-dimensional graded vector space, 
we instead have the space of sections T of a graded vector bundle E over a smooth manifold 
M. We leave the regularity of the elements of £ unspecified (i.e. whether they are smooth 
or distributional), it usually being clear from the context. The dual space £* of the smooth 
elements of £ consists of distributions on M valued in the dual bundle of E. Multilinear 
functionals on £ are elements of T* (g) • • • (g) T*, where the tensor product is completed^^. 

Given v G £, one defines dy : above. Likewise for K G which 

we regard as an integral kernel of an operator (defined with respect to a pairing on £, see 
Section 3.3), we can define Bk as above. If K is not smooth, dx may be ill-defined when 
evaluated on a multilinear functional. In particular, when K arises as the integral kernel of 
a differential operator, one can obtain divergent Feynman integrals from the contractions 
dx applied to local functionals. An element of multilinear map (^J*)®” is said to be local if 
is given by the integral of a polydifferential function of £ over M. 


Appendix B. Wick’s Theorem 


Wick’s Theorem comes in two cases: bosonic and fermionic. The first case gives us a 
combinatorial formula for the integration of monomials against Gaussian measures. We 
assume the reader is familiar with this result and only record it here for notational pur¬ 
poses. Consider with the standard inner product (•, •) and let A = Aij be a symmetric 
nondegenerate dxd matrix. It determines a bilinear form (•, A-) and a normalized Gaussian 

measure dfiA = 

Lemma B.l. (Bosonic Wick’s Theorem) Consider the monomial f{x) = Then 


j dnAf{x) 


1 

2 ”n! 


y 4 L(l)L( 2 ) 

0'^S2n 


. . . ^^cr{2n —1 )^(t(27t) 


(B.l) 


We have x M 2 ) = ® (7°°(M 2 ), where the right-hand side is completed in the sense 

of nuclear Frechet spaces. This readily extends to the case of smooth sections of a vector bundle, so that 

T{EimE2) =V{EA®T{E2). 
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Here denotes the inverse matrix of Aij. 

As is well-known, the sum on the right-hand-side has an elegant description in terms of 
Feynman diagrams. For further details, see e.g. [15, 34]. 

Next, consider the fermionic case of Wick’s Theorem. For this we have to introduce the 
notion of integration over fermionic (odd) variables. Let V be an odd vector space spanned 
by ^ 1 ,... Then we can define the partial integration operator f df,i by 

J dfi = %, 

where 9^. is defined as in (A.l). Thus, 

J dUi = 1 

and more generally 

J di.uiO = fiO 
J dUiO = 0 . 

if /(^) does not depend on We write f df^d ... as shorthand for f d^d ''' / d(,i- Thus, 

J dfd - ■ ■ dfifi.. .f,d = 1. 

Suppose we have an even number of Grassman variables d = 2m. Given a nondegenerate 
skew-symmetric matrix we get the bilinear expression which we abbreviate 

as (^, Aff)/2. Letting d/r = ... dfd-, we have the following: 

Lemma B.2. We have 

J ^ pfaf(A). 

Let 

diXA = Pfaf(A)-id;ue-(^’^«)/2 

denote the normalized “Gaussian measure” on odd variables. We can integrate (polynomial) 
functions in the against this density. We have the identity 

- ^,/(e) + A*"%/(0) 

Letting /(^) = we conclude that 

j d^iA^iCj = A'T (B.2) 

Iteration of this yields the following formula: 

Lemma B.3. (Fermionic Wick’s Theorem) Consider the monomial /(^) = fii ■ ■ ■ f,i 2 „ ■ 
Then 


d^AfiCj = ^ (_1)'^^G(1)V(2) . . . ^G(2n-l)G(2n)_ 


(B.3) 
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It is often the case that the set of odd variables is partitioned into two separate sets, 
and their corresponding “conjugate” variables cjj,..., Moreover, we are 
given the bilinear expression cj* BijUj with Bij is an arbitrary matrix, which we abbreviate 
by u:*Buj. One can think of this as letting and ^m+k = for 1 < /c < m and 



(B.4) 


in the above. Let dfi = dcoldLOi ■ ■ ■ duj'^dojm- 

Lemma B.4. 

J d^e~‘^ = det B. 

Since 

Pfaf(A) = (-l)("2 ")det(B) 

/2m\ 

d^i - ■ ■ d^2m = (- 1 )^ 2 /duldui • • • dco^dum, 

the formulas in Lemmas B.2 and B.4 agree. 

We suppose B^j is invertible. Let dfiB = From 

we have 

J diiBOJiijj] = , 

or in other words 


J djiB^i^j = — P*-^ . 


(B.5) 


One can attribute this minus sign to the minus sign that occurs in 



(B.6) 


for A of the form (B.4). 

Remark B.5. The above considerations explain the convention that fermionic loops in 
quantum field theoretic computations are weighted with a minus sign. This is a basis 
dependent statement, however. One should regard (B.3) as fundamental, with (B.5) a 
consequence of a particular (albeit common) parametrization. 

One can unify the bosonic and fermionic cases of Wick’s Theorem using the framework 
of graded vector spaces. Let A{-, ■) be a symmetric pairing on a graded vector space V, i.e., 
an element of Sym^(P*), which is nondegenerate. Such a pairing determines a dual pairing 
P on V* which is an element of Sym^(y). We call P the propagator. 

If we pick a basis Vi (of homogeneous degree elements) for V, then if Aij = A(vi,Vj), we 
have P = A^^ViVj. Define 

dtiA = CA n 

with CA chosen so that so that f d^A = 1- We have the following; 
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Lemma B.6. (Wick’s Theorem, unified version) For f a (polynomial) funetion on V, 

I dfiAfix) = (e^^/)(0). (B.7) 

Here, dp : Sym"'(H*) —)■ is the Wick contraction operator, defined as in 

Section A.l, and 

_ V jdpT 

ho 

is the sum over all Wick contractions weighted with the appropriate symmetry factor. 
The right-hand side of (B.7) is evaluated at zero so that the maximal number of Wick 
contractions have been made (in the Feynman-diagrammatic picture, we only sum over 
vacuum diagrams, i.e. those without external legs). 
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